Light propagation in non-linear electrodynamics 
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Abstract 

Working on the approximation of low frequency, we present the light cone conditions for a class of 
theories constructed with the two gauge invariants of the Maxwell field without making use of average 
over polarization states. Different polarization states are thus identified describing birefringence 
phenomena. We make an application of the formalism to the case of Euler-Heisenberg effective 
' Lagrangian and well know results are obtained. 

>V PACS numbers: 42.25.L, 41. 20. J 
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CO ■ I. INTRODUCTION 

It is now a well established result that the velocity of the electromagnetic waves has its value dependent on the 
. vacuum polarization states. Indeed, such polarization effects appear when a strong field (electrodynamical critical 
Q\ ' field: E ar = B cr = m 2 c 2 jeh « 1.3 x 10 ls V/m as 4.4 x 10 13 G) is produced in some region of space. The most important 
consequence of this fact consists in the birefringence effect: the velocity of wave propagation depending on the wave 
polarization. A experimental method to detect the vacuum birefringence induced by a magnetic field was proposed 
in 1979 by E. Iacoppini and E. Zavattini In the experimental context, it is worth to mention the work of D. 
Bakalov et al 0, where optical techniques is used to detect birefringence in the presence of a strong magnetic field 
(PVLAS experiment). The theoretical description of nonlinear effects on light propagation was studied long before 
by Z. Bialynicka-Birula and I. Bialynicki-Birula ||, where was calculated the probability of the photon splitting in 
an external electromagnetic field. The same problem was extensively studied by S.L. Adler j|. Other beautiful new 
I , results on vacuum polarization phenomena in nontrivial vacua, including curved spacetimes, can be found in the works 
of J.I. Latorre, P. Pascual and R. Tarrach j§ , I.T. Drumond and S.J. Hathrell §|, G.M. Shore @ and others. See 
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for instance the works of K. Scharnhorst and G. Barton where the problem of photon propagation between 
parallel mirrors is worked out. 

Recently, W. Dittrich and H. Gies p0| , within the approach of the geometrical optics, derived the light cone 
conditions for a class of homogeneous nontrivial QED vacua using the rule of average over polarization states. They 
generalized some results previously obtained by Latorre and others, in particular the so called "unified formula" ||. 
Indeed, such unified formula was identified by Latorre et al for several modified QED vacua and proved by Dittrich 
and Gies for certain cases. The use of the above mentioned average procedure excludes from their formalism the 
possibility to analyse the important phenomena of birefringence. 

In this paper we deal with a class of Lagrangians depending on the two Lorentz and gauge invariants of the Maxwell 
field 

L = L{F^F flv ,F^F flv ) (1) 

and working on the approximation of soft photons (the wavelength of propagating wave is large compared to the 
Compton wavelength) , we present the light cone conditions for local theories of gauge invariant spin-one fields without 
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making use of average over polarization states. Recent results by Dittrich and Gies are thus generalized in our 
approach, whose main contribution relies on the study of birefringence effects in a unified formalism. The polarization 
problem is worked out and the dispersion law is obtained, showing that there are a different polarization mode 
associated to each velocity of wave propagation. Finally, we apply the formalism to Euler-Heisenberg Lagrangian and 
some known results are derived in the context of the present formalism. We set the units c = 1 = h. 



II. NON-LINEAR SPIN-ONE THEORIES 



Instead of calculating light cone conditions for each particular theory, we make use here of a general formalism, 
applicable to any Lagrangian based local theory describing gauge invariant spin-one fields that can be constructed 
with the two invariants of the Maxwell field. We denote the electromagnetic field strength by the anti-symmetric 
2-rank tensor , and its dual is defined as 

F a f3 = \v a!3 ° T F aT . (2) 

Let us set the only two local and gauge invariant scalar fields F and G associated with by 

F = F^F^ (3a) 

G = F^F^. (3b) 

In order to achieve more simplicity, we work in Minkowski spacetime employing a Cartesian coordinate system. 
Thus, the background metric will be represented by r/^, which is defined by diag(+l, —1, —1, —1). We defined the 
completely anti-symmetric tensor r^P^ (^ 0123 = 1), and set the notation Lx = BL/dX , where the variable X stands 
for any monomial on the field invariants. The gauge invariant density of Lagrangian of electrodynamics is an arbitrary 
function of F and G: 



L = L(F,G). (4) 
From the minimal action principle we get the equation of motion 

^L F F" V + L G F^ ,„ = (5) 

where a comma denotes partial derivatives with respect to the Cartesian coordinates. Using relations F v — 2F a i 3 F a p ll 
and G M = 2F al3 F " a f3,v in equation (||) we obtain: 

2N^F a0 , u + L F F^, v = O (6) 

where we introduced the 4-rank tensor Nv a P through 

N" uaf} = L FF F^F aP + L GG F^F al} + L FG (f^F " 3 + F^F af ^ . (7) 

Additionally, the field strength F^ must satisfy the Bianchi identity. 

Let us now turn our attention to the expressions which represent the light cone conditions for such general non- linear 
spin-one theory. 



III. WAVE PROPAGATION 



In this section we analyse the propagation of linear shock-waves associated with the discontinuities of the field in 
the limit of geometrical optics O]. Let us consider a surface of discontinuity S defined by 

z(x") = 0. (8) 

Whenever S is a global surface, it divides the spacetime in two distinct regions U~ and U + (z < and z > 0, 
respectively). Given an arbitrary function of the coordinates, /(a; M ), we define its discontinuity on S as 
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[/(x«)] s = {p l ± im p [/(P+) - f(P~)] (9) 



where P + , P and P belong to U + , U and £ respectively. Applying the conditions |11| of discontinuities for the 
tensor field P M „ and its derivatives we set 

[P Q/3 ] S = (10a) 
[FapAn = f<*pk\ (10b) 

where f a p represents the discontinuities of the field on the surface £ and k\ is the wave propagation 4-vector. The 
discontinuity of the Bianchi identity yields 

faflk\ + ff3\k a + fxak/3 = 0. (11) 

* 

In order to obtain scalar relations, we consider the product of equation ( pd| ) with F a,3 k x and with F a @k x , which leads 
to 

(k x k x = -2F a ^f px k x k a (12a) 
Ck X k x = -2F a/3 f px k x k a (12b) 
where we have introduced the scalar quantities C and C* as 

C = F^f aP (13a) 
C = *F al} .f a p. (13b) 
In the same way, taking the discontinuities of the field equation (||), we get 

fpxk x = -^-N ( r p f vp K. (14) 
hp 

Introducing relation ( O ) in equations ( |l2| ) and making use of the useful identities between the tensor field F M „ and 
its dual: 

F m F a v = —Grifiv (15a) 



4 

1 

r 



F m F» a - F^F» a = -F8S (15b) 



we obtain 



(k 2 = -^-F^F\k v k T (L FF ( + L FG C) - ^-k 2 (L FG ( + L GG C) (16a) 

A /~1 Of? 

C * fc 2 = F ^F\k u k T (L FG ( + LggC) ~ -^k 2 (L FF ( + L FG C) + ^k 2 (L FG ( + L GG C) (16b) 
L/ F L/ F L/ F 

where k 2 = rj^k^k^. 

In order to find the equation that represents the propagation of the field discontinuities, we seek for a master 
relation which should be independent of the quantities f a p (that is, independent of £ and £*). There is a si mple way 
t o ac hieve such goal. We firstly isolate the common term P T pk v k T which appears in both equations (16a) and 
(|16b|) . The difference of these equations can thus be written as 



(k 2 Ck 2 2Fk 2 Gk 2 (L ff C + L F gC L fg { + L G gC 



L FF ( + L FG (* L F gC + LggC* L F hp \Lp G ( + L GG (* L FF C, + L FG (* 

Assuming k 2 ^ 0, we obtain an algebraic linear relation 



(17) 



1 Note that k 2 = represents the standard propagation which occurs for the linear theory L = —F/4. We are interested here 
in the study of the possible deviations from this simple case. 
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(18) 



between £ and , where we denned 
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Lfg + 2 7 — LfqLqq - —L 2 FG + — L 2 GG 

Lip Lip Lip 



F 



F 
'Lp 



Lgg — Lff + 2- — LfpLgg + 2- — L FG — 2- — LffLfg + 2- — LfgLgg 



G_ 



(3_ 

Lp 



Lfg + "2- — LffLfg — 7 — L 2 FF + - — L FG - 
Lp hp Lip 



Solving the quadratic equation (|lj) for £* we obtain 



with 
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-n 2 ± yjvi 2 . - 
2^ 



(19a) 
(19b) 
(19c) 

(20) 

(21) 



Using this solution into equation ( |16a| ) and assuming £ ^ it results the following light cone conditions for spin-one 
fields: 



G 

1 + — {L FG + n ± L GG ) 
Lp 



k 2 - — {Lff + ^±L FG ) F^F\k u k T = 0. 



In a similar way for equation ( 16b ) we obtain 



2F G 
—— {Lfg + Q±Lgg) + 7 — {Lff + ^1±Lfg) 

L p Li F 



k 2 - — {L FG + Sl ± L GG ) F^F\k v k T = 0. 



(22) 



(23) 



After some algebra, it can be shown that equation (g3|) is identical to (p2|). 

The light cone conditions following from such procedure lead to two possible paths of propagation, according with 
the double solutions £1 ± . These two conditions are related to distinct polarizations modes, as we will see in the next 
section, indicating the possibility of birefringence. This effect depends upon the particular theory we shall consider. 

In a more appealing form, we can present expression (E2h as 



k = 4 



Lff + SI±Lfg 



Lp + G {Lp G + Vt ± L GG ) 



F x ^F v x Kk v . 



Since k 2 = lo 2 — |fc| 2 , the phase velocity w /|A;| = v of the propagating light is found to be 



v* = 1 - 4 



jFF 



n + L 



±L,p G 



Lf + G {Lfc + £1±Lgg) 



F^F\k^k v . 



(24) 



(25) 



This equation indicates that the familiar case k 2 = 0, which occurs for linear electrodynamics L = —F/A, is also 
possible for more complicated situations, as for those theories satisfying the relation Lff + Q±Lfg = 0. Solutions of 
this equation brings up the form of particular Lagrangians for which k 2 = 0, despite it was previously assumed that 
k 2 5^ 0. A simple solution of this equation can be obtained by setting L = — F/A + f{G), where f{G) is an arbitrary 
function of the invariant G. Another example where k 2 — occurs consists in the nonlinear Lagrangian of the N. 
Born and L. Infold |12[ . We are not interested here on the analysis of such theories. 

The light cone conditions (E4h can also be expressed in terms of the energy-momentum tensor of the non-linear field 



7) t „ = -AL F F^ F av — {L — GL G ) r,^ 



(26) 
(27) 



where we defined the quantities 
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Q ± = — (L FF + n ± L FG ) (2g) 

Lp + GLp (Lpc + Q±Lgg) + {Lff + Q±Lfg) (L — GLq) 

Thus, in terms of the energy momentum tensor, the phase velocity for each propagation mode, corresponding to the 
solutions of fi ± , are given by 

4 = 1- Q ± T^ ntl n v (29) 

where we introduced the quantity = k^/\k\ specifying the direction of wave propagation. 

In the literature one usually makes use of an average over polarization states, which can directly be obtained from 
our formalism in the form v = (v + + u_)/2. This represents the velocity of propagation for the average mode. 

The results we had obtained in this section are applied for spin-one theories which is set by the Lagrangian function 
defined by (H), in the approximation of soft photons. The light cone conditions, here presented in terms of the field 
strength by (pif) or else in terms of the energy- momentum tensor by (^), are useful in a variety of situations, and 
particularly for the study of birefringence phenomena. Let us now analyse in more details the problem of polarization. 



IV. POLARIZATION 



The most general decomposition for a skew-symmetric tensor is f a p — {A a Bf3 — ApB a ) + (C a Dp — CpDa), where 
t he v ectors A a , B a , C a , D a are arbitrary. For the case where f a p is the wave propagation tensor given by equation 
(10t), for which equation (O) applies, it follows that the above decomposition simplifies to 



f a p = a(e a kf3 - epk a ) 



(30) 



where a is the strength of the wavelet and ep represents the polarization vector. 

The field equations impose restrictions on the possible states determined by such vector. Introducing (|3^) in ( |l4| ) 
and assuming a/0we have 



k 2 e" = — — (N^Kk^e". 



(31) 



The k 2 = case includes the linear propagation regime, where the polarization modes are well known. In this case 
equation © states that uf3 k a k fi e u = 0. The k 2 j= case can be treated defining a symmetric tensor Z M „ by 

4 



Lpk z 

With it we can write ( |3l| ) as an eigenvector equation 



N» a „ k a k . 



0. 



(32) 



(33) 



The solutions of equation ( p3|) (eigenvectors of Z M „) represent the dynamically allowed polarization modes. The 
definition of N^ 01 ^, from equation (0), leads us to conclude that the tensor structure of Z^ v can be completely 
determined by the electromagnetic tensor, its dual and the wave vector Hence, the general solution for the 
eigenvector problem can be achieved by expanding e M as a linear combination of the following linearly independent 
vectors: 



F^k v ee o", F^K = a M , F^ a F av k v = b", k". 
Thus, the polarization vector takes the form 



(34) 
(35) 



Introducing ( |35| ) in (33) and taking the products, results 



— k z + Lppa + LpGd^a^ 



a 



L FG a 2 + Lcca^a, 



■0 



Lppa^a^ + LpQa^a^ 



+ -f[0]k u + 6 



^k 2 
4 



b" = 0. 



(36) 
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In order to satisfy the above equation, the coefficients of each independent vector must be null, resulting in 



F 9 * 

— k + L FF a + Lpoa ra^ 



+ (3 



L FG a + Laaa^a, 





7 = arbitrary 
5 = 0. 



(37a) 

(37b) 

(37c) 
(37d) 



Therefore, from 
P) yields 



]), 7/^ does not contribute to f a p, and we shall not consider it in any further. Using the relations 



* u * 2 

& (Z^ = a 



iFk 2 . 



(38a) 
(38b) 



Substituting these results in equations (^) and solving the system for k 2 we obtain the dispersion laws j|[l4j for both 
polarization modes, which are described by the vectors e±. Indeed, it can be shown that such dispersion relations are 
recovered from the light cone conditions (^4|) for all cases known in the literature, ensuring our previous statement 
concerning the relation between f2 ± with the two polarization states . 

V. APPLICATION TO EULER-HEISENBERG LAGRANGIAN 

In this section we apply the previous results to retrieve, from our formalism, a particular case of birefringence 
presented in the literature || . The most familiar non- linear case of electrodynamics is given by Euler-Heisenberg 
effective action of QED^|, 



S= dx 



- l -F+ l - [ F 2 + -G 2 
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where dx stands for the volume measure of the spacctime, and the quantum parameter fi is defined by 

2a 2 



45to 4 



(39) 



(40) 



The light cone condition for each propagation mode can be directly obtained from expression ( |24| ) as m 

k 2 = -8fxF a ^F/k a k (44a) 
k 2 = -UixF at *F/k a k . (44b) 

From equation (p8|) , in the required order of approximation we obtain 

(42a) 
(42b) 



0+ = Sfi 
Q_ = 14/i. 

Thus, in terms of the energy-momentum tensor, 

k 2 = -8/iT Q/3 fc a fc /3 
k 2 = -UnT al3 k a kf3. 

The phase velocities ( p5|) corresponding to ( |44| ) reduce to 



(43a) 
(43b) 



The limit of weak-field from one loop approximated QED action. 



G 



v+ = 1 - ^F a ^F^n a n p (44a) 
v- = 1 - 7fj,F a,1 F/n a np. (44b) 

These velocities correspond to two polarization states, showing that birefringence effect occurs in Euler-Heisenbcrg 
model, as it is well known in the literature. From the polarization sum rule we obtain 

v 2 = 1 - ll/iF^F^n^. (45) 

This result was presented in the paper of G. M. Shore Q, and more recently by W. Dittrich and H. Gies JlO[ . 



VI. CONCLUSION 



In this work we derived light cone conditions for a class of local and gauge invariant spin-one field theory constructed 
with the two invariants of the Maxwell field in the approximation of low frequency without making use of the average 
over polarization modes. In our formalism we took the different polarization states into account explicitly and one 
achieved a generalization of the wave propagation formula, which can be applied to the analysis of birefringence 
effects in a unified way. In order to illustrate its applications, we exhibited how to obtain some results of quantum 
electrodynamics concerning the wave velocity dependence on polarization states. 

An interesting continuation of this work would be the analysis of the wave propagation equations as an effective 
modification of Minkowskian geometry for each polarization direction, due to non-linear effects. Investigations on this 
topic had already been performed (see reference (l5"| , and more recently [[l6||). The use of the formalism presented 
here to the case of wave propagation in non-linear material media, also deserves further investigations. 
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